Symbiotic gap and semi-gap solitons in Bose-Einstein condensates 
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Using the variational approximation (VA) and numerical simulations, we study one-dimensional 
gap solitons in a binary Bose-Einstein condensate trapped in an optical-lattice potential. We consider 
the case of inter-species repulsion, while the intra-species interaction may be either repulsive or 
attractive. Several types of gap solitons are found: symmetric or asymmetric; unsplit or split, if 
centers of the components coincide or separate; intra-gap (with both chemical potentials falling into 
a single bandgap) or inter-gap, otherwise. In the case of the intra-species attraction, a smooth 
transition takes place between solitons in the semi-infinite gap, the ones in the first finite bandgap, 
and semi-gap solitons (with one component in a bandgap and the other in the semi-infinite gap). 

PACS numbers: 03.75.Ss,03.75.Lm,05.45.Yv 



I. INTRODUCTION 

One of the milestones in studies of Bose-Einstein con- 
densates (BECs) was the creation of bright solitons in 7 Li 
and 85 Rb in "cigar-shaped" traps [l|, with the atomic 
scattering length made negative (which corresponds to 
the attraction between atoms) by means of the Feshbach- 
resonance (FR) technique 0. Normally, BEC features 
repulsion among atoms. In that case, it was predicted 
that an optical-lattice (OL) potential may support gap 
solitons (GSs) [|| , whose chemical potential falls in finite 
bandgaps of the OL-induced spectrum. Although GSs, 
unlike ordinary solitons in self-attractive BEC, cannot 
realize the ground state of the condensate, it was demon- 
strated that they may easily be stable against small per- 
turbations 0. A GS in Rb was experimentally created 
in a cigar-shaped trap combined with an OL potential, 
pushing the BEC into the appropriate bandgap by accel- 
eration Q . Other possibilities for the creation of GSs are 
offered by phase imprinting or squeezing the system 
into a small region by a tight longitudinal parabolic trap, 
which is subsequently relaxed Q- 

BEC mixture of two hyperfine states of the same atom 
are also available to the experiment [|[. The sign and 
strength of the inter-species interaction may also be con- 
trolled by means of the FR Q , hence one may consider a 
binary condensate with intra-species repulsion combined 
with attraction between the species. It was proposed to 
use this setting for the creation of symbiotic solitons [Icl ] , 
in which the attraction overcomes the intrinsic repulsion. 

In this work, we aim to study compact (tightly bound 
pd| ) symbiotic gap solitons in a binary BEC, which are 
trapped, essentially, in a single cell of the underlying OL 
potential. Unlike the situation dealt with in Rcfs. [l(|, 
we consider the case of inter-species repulsion, while the 
intra-species interactions may be repulsive or attractive. 
In Ref. [ll| it was already demonstrated that the addi- 
tion of intra-species repulsion expands the stability re- 
gion of symbiotic GSs supported primarily by the inter- 
species repulsion. The case of attraction between two 
self- repulsive species was recently considered in Ref. fl2| , 



where it was shown that the attraction leads to a counter- 
intuitive result - splitting between GSs formed in each 
species. This effect can be explained by a negative effec- 
tive mass, which is a characteristic feature of the GS [|[. 
Indeed, considering the interaction of two GSs belonging 
to different species, one may expect that the interplay 
of the attractive interaction with the negative mass will 
split the GS pair. 

Using variational [l3T ] and numerical methods, we here 
construct families of stable GSs of two kinds: unsplit 
(fully overlapping) and split (separated). The split- 
ting border is predicted by the variational approximation 
(VA) in an almost exact form. In terms of chemical po- 
tentials of the two components, the solitons may be of 
intra- and inter-gap types [Tl| , with the two components 
sitting, respectively, in the same gap or different gaps. In 
particular, the states with one component residing in the 
semi-infinite gap (which is possible in the case of intra- 
species attraction) will be called semi-gap solitons. 

The paper is organized as follows. The formulation 
of the system and analytical results, obtained by the 
variational method [l3j . are given in Sec. II. Numeri- 
cal findings are reported in Sec. Ill, including maps of 
GS families in appropriate parameter planes. Section IV 
summarizes the work. 



II. ANALYTICAL CONSIDERATIONS 

We consider a binary BEC loaded into a cigar-shaped 
trap combined with an OL potential acting in the axial 
direction. Starting with the system of coupled 3D Gross- 
Pitaevskii equations (GPEs) for wave functions of the 
two components, 4>\ and <f>2, one can reduce them to ID 
equations In the scaled form, they are [l2j 

*(^,a) t = -(1/2)^1,2)^+ <?|<M 2 </>i,2 

+ .912 |02,1 | 2 <?H,2 ~ V COS (2x) 01,2, (1) 

where the OL period is fixed to be n, and the wave func- 
tions are normalized to numbers of atoms in the two 
species, f^°° \<j>\ t 2(x) \ 2 dx — Ni t 2- In Eq. |T]), time, 



2 



the OL strength, and nonlincarity coefficients are related 
to their counterparts measured in physical units as fol- 
lows: t = (ir/L) (h/m)t phys ,V = (L/nh) m(Vb) p hys, 
{9t9i2} = (2Lmu}_i_/irh) {0,012}, where m is the atomic 
mass, L the OL period, a and 012 scattering lengths 
accounting for collisions between atoms belonging to 
the same or different species, and lo± the transverse- 
confinement frequency. As said above, we assume re- 
pulsive inter-species interactions, with gi 2 > 0, while the 
intra-species nonlinearity may be both repulsive (g > 0) 
and attractive (g < 0). 

While the model assumes equal intra-species scatter- 
ing lengths, they are, in general, different for two hy- 
pcrfine states [8|. Therefore, using a FR, one cannot 
modify both intra-species nonlinearities to keep exactly 
equal values of coefficient g in equations for both com- 
ponents [cf. Eq. ([Tjl], running from negative to positive 
values (hence, strictly speaking, different cases consid- 
ered in this work cannot be realized in a single mixture, 
but should be rather considered as a collection of sit- 
uations occurring in different mixtures). However, we 
will consider asymmetric configurations, with N\ ^ N 2 , 
which give rise to a much stronger difference in the effec- 
tive interaction strengths in the two components than a 
small difference in their intrinsic scattering lengths. 

Stationary solutions to Eqs. |T]) are looked for in the 
usual form, (j)i^ 2 (x 1 t) = exp (— i/ii^i) u\^ 2 {x), with chem- 
ical potentials /ii.2 and functions ui^(x) obeying 



to Ni,2- Thus, the remaining equations, dL/dwi^ 
<9-L/<9Ki,2 = 0, take the form 
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Using Eqs. (J5j) and j6]) we can predict borders between 
intra-gap and inter-gap soliton families of different types. 
To this end, we take /j,\ i2 from Eqs. ^ and, referring 
to the spectrum of the linearized equation (fT]), identify 
curves in plane (iVi, N 2 ) which correspond to boundaries 
between different gaps in the two components. 

Variational approximation for split solitons: Two- 
component solitons different from those considered above 
feature splitting between the two components. An is- 
sue of obvious interest is to predict the splitting thresh- 
old by means of the VA, for the symmetric case, with 
Ni = N 2 = N. For this purpose, we use the following 
ansatz, 



(split) 
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A»i,2Wi,2+«i 2/2- gu\ 2 - g\ 2 u\ cos (2x) u 1)2 = 0, 

(2) 

with /_ u 2 2 (x)dx = Ni )2 . In the GS solutions con- 
structed below, fi\ and fi 2 belong to the first two finite 
bandgaps and/or the semi- infinite gap in the spectrum 
induced by potential — Vq cos (2x). 

Variational approximation for unsplit solitons: Equa- 
tion © can be derived from Lagrangian 



with infinitesimal splitting parameter 6, the objective 
being to find a point at which a solution with 6^0 
emerges. At small b, the two components of expression 
(J?]) feature maxima shifted to x = ±b/a + 0(b 2 ), and 
up to order b 2 , it satisfies the normalization conditions, 
u i 2 {x)dx = N . Unlike b, constant C, to be defined 
below, is not a variational parameter. 

The substitution of ansatz |(7j) in Lagrangian © yields, 
at orders b° and b 2 . 



L 



Hiul + [i 2 u\ - - 
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+V cos(2x)(u 2 +ul)- -g (uf + uf) 
-gi 2 u 2 ul] dx - H1N1 - n 2 N 2 . 



(3) 



To predict solitons with a compact symmetric profile, 
which corresponds to numerical results displayed below, 
we adopt the Gaussian ansatz [l3| . 



v 
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where variational parameters are widths u>i.2, reduced 
norms Hi^, and ^1,2- The substitution of the ansatz 
in Eq. ([3]) yields an effective Lagrangian, L = 
L (^i,2) Wi,2, Mi, 2)- Then, the first pair of the variational 
equations, dL/dfii, 2 = 0, gives Hi. 2 = 1, which is substi- 
tuted below, after performing the variation with respect 
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At order 6° (i.e., for the unsplit soliton), variational equa- 
tion dL/dw = reduces to Eq. © with N ± = N 2 = N 
and w\ = w 2 = w: 



1 + (g + g ^ iVw = W w*e- 2 



'2-k 



(9) 



At order b 2 , equation <9L/d(o 2 ) = yields the splitting 
condition. 



gNw\ C~2g 12 Nw 



CV w 4 e' 1 "- = 0. 

(10) 

Obviously, the splitting should not occur if g\ 2 = 0, i.e., 
Eq. (TIT))) must only yield the trivial solution, w = 0, 
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FIG. 1: (Color online) Examples of unsplit and split sym- 
metric solitons, with Ni = N% = 1000, trapped in potential 
— Vbcos(2x). Here and in all other figures, Vo = 5. For the 
unsplit soliton, the variational profile is included too. 



in this case. This condition selects the value of C which 
was arbitrary hitherto: C = —2, hence Eq. (p~0|) takes the 
form Ngi2 = 2\f2TtVow 3 e~ w . Combining this with Eq. 
©, we obtain w = ^/2tt/ [iV(<7i2 — g)], and a prediction 
for TV at the splitting point: 



8tt 2 Vo 



512 (512 - g)' 



exp 



2tt 



(312 - 9) 



(11) 



III. NUMERICAL RESULTS 

Symmetric solitons: Equation |T]) was discrctized us- 
ing the Crank-Nicholson scheme and solved numerically 
in real time, until the solution would converge to a sta- 
tionary soliton. This way of generating the solitons 
guarantees their stability. In Fig. [TJ we present typi- 
cal profiles of split and unsplit symmetric solitons, with 
Ni = N-2- Due to the symmetry, these solitons are al- 
ways of the intra-gap type (in Fig. [H they belong to the 
second bandgap; in the semi-infinite and first gaps, the 
shape of the solitons arc quite similar). 

In all cases, the difference between the variational and 
numerical shapes of the unsplit solitons is extremely 
small. The present solitons are essentially confined to 
a single cell of the OL potential. They change the shape 
and develop undulating tails, which are often considered 
as a characteristic feature of GSs, when [i is taken very 
close to an edge of the bandgap (Fig. [1] demonstrates 
that, even in a well-pronounced split state, peaks of both 
components stay in a common cell). It is also observed 
that, as might be expected, the increase of the intra- 
species nonlinearity coefficient, g, pushes the solitons to 
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FIG. 2: (Color online) The family of the two-component sym- 
metric solitons (Ni — N2 = N), mapped into the plane of the 
interaction strengths, g^N and gN. The plane is divided 
into regions corresponding to the semi-infinite gap and the fi- 
nite first and second bandgaps. They are separated by narrow 
stripes representing the Bloch bands. The border between the 
unsplit and split solitons is shown as found from the numerical 
data, and as predicted by Eq. (TiTj) . 



higher bandgaps, while the increase of gi2 tends to split 
the two components of the soliton. In addition to the 
compact GSs presented here, there may also exist loosely 
bound ones, that extend over several OL [la) . 

In Fig. O the entire family of the symmetric solitons 
is displayed in the parameter plane of the strength of the 
intra- and inter-species interactions, (gN, gnN). The VA 
prediction for border between the unsplit and split soli- 
tons, given by Eq. (jlip , provides a remarkably accurate 
fit to the numerical findings. 

Dependences N(fi) for families of symmetric solitons 
are plotted in Fig. [3J It is known that a necessary sta- 
bility condition for solitons populating the semi-infinite 
gap is given by the Vakhitov-Kolokolov (VK) criterion, 
dN/dfi < [la ], while stable solitons in finite bandgaps 
have dN/d/j, > 0, disobeying this criterion In the 

present case, Fig. [3] shows the same generic feature (the 
semi-infinite gap contains solitons only for g < 0, i.e., 
in the case of the self- attraction). A noteworthy feature, 
viz., a turning point in dependence N(/i), is exhibited, for 
g = —0.001, by the solution branch which passes from 
the semi-infinite gap into the first finite bandgap, and 
also by the branch corresponding to g = —0.0005. Con- 
sequently, two different stable solitons can be found in 
the corresponding interval of /i. The solitons belonging 
to the branches with g = 0.0025, g = 0.001 and .9 = 
in Fig. [3Jare unsplit, and they are accurately predicted 
by the VA. Accordingly, the curves for these branches, 
as obtained from the VA and from the numerical data, 
are virtually identical. On the other hand, all solitons be- 
longing to the branch with g = —0.0025 exhibit splitting. 
As concerns the bending branches, their parts below the 
turning point are formed by unsplit solitons (which are 
accurately approximated by the VA), while above the 
turning point the family continues in the split form. Ac- 
cordingly, the turning point on each bending branch be- 
longs to the splitting border for the symmetric solitons, 



FIG. 3: (Color online) The number of atoms in the symmetric 
soliton, Ni — N2 = N, versus the common chemical potential 
of both components, at several values of g for gi2 = 0.002. 
Vertical stripes are the Bloch bands between the gaps (the 
solution branch with g = —0.001 suffers a discontinuity when 
it hits the band separating the semi-infinite and first finite 
gaps). 
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FIG. 4: (Color online) Typical profiles of unsplit and split 
asymmetric (Ni 7^ N2) solitons. The examples represent soli- 
tons of inter-gap types, as indicated in the panels. For the 
unsplit soliton, the profiles predicted by the VA are shown 
too. 



cf. Fig. m 

Asymmetric solitons: Typical examples of solitons 
with Ni ^ N2 are displayed in Fig. [4j Similar to their 
symmetric counterparts, cf. Fig. [H they feature both 
unsplit and split shapes (the former ones are well ap- 
proximated by the VA), which are again confined to a 
single cell of the OL potential. 

The entire family of asymmetric and symmetric GSs 
is mapped in the (Ni, N2) plane, at fixed values of the 
interaction coefficients (c/12 and <?), in Fig. [5] In these 
diagrams, the border between intra-gap solitons of differ- 
ent types shrink to a point belonging to the diagonal line 
(Ni = N2), which corresponds to symmetric solitons that 
account for direct transitions between different types of 
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FIG. 5: (Color online) Families of asymmetric and symmetric 
solitons mapped into the plane of atom numbers iVi and N3, 
for different gi2 and g. The plane is divided into regions 
populated by solitons of six different types (three intra-gap 
and three inter-gap varieties, symbols and 1, 2 standing for 
the semi-infinite and two lowest finite bandgaps, respectively). 
Each panel also shows the numerically found border between 
the unsplit and split solitons, and borders between different 
types of the unsplit ones, as predicted by the VA. 



intra-gap solitons. In Fig. [5] the VA for the unsplit soli- 
tons accurately predicts borders between their different 
varieties. 

If none of the nonlinearities is attractive [Figs. EJa) 
and (b)], no chemical potential may fall in the semi- 
infinite gap. Three types of GSs are possible if both 
nonlinearities arc repulsive [Fig. EJa)]: intra-gap ones, 
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in the two finite bandgaps, and the inter-gap species, 
combining them. If the intra-species nonlinearity exactly 
vanishes [Figs. E|b)], the inter-species repulsion cannot 
push both components into the second finite bandgap, 
which leaves us with two species: intra-gap in the first 
bandgap, and the one mixing the two finite bandgaps. 
The interplay of the attractive intra-species nonlinearity 
with the inter-species repulsion supports two intra-gap 
and two inter-gap types, as seen in Fig. [5jc). Note that 
one of them skips the first bandgap, binding together 
components sitting in the semi-infinite and in second fi- 
nite gaps. A notable feature of the map in Fig. [5jc) is 
the smooth transition from ordinary solitons, with both 
components in the semi-infinite gap, to ones of the semi- 
gap type. 

IV. CONCLUSION 

In this work, we have considered the interplay of the 
repulsion between two species of bosonic atoms with 
intra-species repulsion or attraction in a binary BEC 
mixture loaded into the OL potential. Families of sta- 
ble solitons found in this setting are classified as sym- 
metric/asymmetric, split/unsplit, and intra/inter-gap. 



Three varieties of intra-gap solitons, and another three 
types of inter-gap ones are identified, if the considera- 
tion is limited to the two lowest finite bandgaps of the 
OL-induced spectrum. Varying the atom numbers in the 
two components, Ni 2> we have plotted maps of various 
states. Although different intra- and inter-gap species 
are separated by Bloch bands, transitions between them 
are continuous in the (Ni,N2) plane. In particular, a 
solution branch which connects the solitons (of the split 
type), populating the semi- infinite gap, and unsplit soli- 
tons in the first finite bandgap, features the turning point 
at the border between the two varieties. Other varieties 
revealed by the analysis represent semi-gap solitons, with 
one component belonging to the semi-infinite gap, and 
the other one falling into a finite bandgap. 

A considerable part of the numerical findings reported 
in this work was accurately predicted by variational ap- 
proximation. These include the shape of unsplit solitons 
(both symmetric and asymmetric ones), borders between 
their varieties, and the splitting border for the symmetric 
solitons. 
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